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In most current models of inflation based on a weakly self-coupled scalar matter field minimally 
coupled to gravity, the period of inflation lasts so long that, at the beginning of the inflationary 
period, the physical wavelengths of comoving scales which correspond to the present large-scale 
structure of the Universe were smaller than the Planck length. Thus, the usual computations of the 
spectrum of fluctuations in these models involve extrapolating low energy physics (both in the matter 
and gravitational sector) into regions where this physics is not applicable. In this article we study 
the dependence of the usual predictions of inflation for the spectrum of cosmological fluctuations 
on the hidden assumptions about super-Planck scale physics. We introduce a class of modified 
dispersion relations to mimic possible effects of super-Planck scale physics, and find that, given an 
initial state determined by minimizing the energy density, for dispersions relations introduced by 
Unruh the spectrum is unchanged, whereas for a class of dispersion relations similar to those used 
by Corley and Jacobson (which involve a more radical departure from the usual linear relation) 
important deviations from the usual predictions of inflation can be obtained. Some implications of 
this result for the unification of fundamental physics and early Universe cosmology are discussed. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The inflationary Universe scenario Q is the first theory 
of the very early Universe to provide a mechanism Q for 
the production of density fluctuations on scales of cosmo- 
logical interest based on causal physics (see also Rcf. || 
for initial ideas). The key point is that during the period 
of inflation fixed comoving scales are stretched exponen- 
tially compared to the Hubble radius. Thus, the wave- 
lengths corresponding to the present large-scale structure 
in the Universe and to the measured Cosmic Microwave 
Background (CMB) anisotropies were equal to the Hub- 
ble radius about 50 Hubble expansion times before the 
end of inflation. This gives rise to the possibility that 
causal physics acting before that time can generate fluc- 
tuations on these scales while they are of sub-Hubblc 
length. 

Most current models of inflation are based on weakly 
self-coupled scalar matter fields minimally coupled to 
gravity. In this context, quantum vacuum fluctuations 
provide (2j a causal mechanism for generating fluctua- 
tions. In fact, the coupled linear metric and matter fluc- 
tuations can be quantized in a unified manner jj]. The 
problem reduces to the quantization of a free scalar field 
with a time-dependent mass (see e.g. Ref. ||] for a com- 
prehensive review). An initial vacuum state thus under- 
goes squeezing during inflation, and this leads to the gen- 
eration of fluctuations. According to the standard calcu- 



lations |2|,|6|-|l0(1 , the predicted spectrum is scale- invariant 
(modulo a mild deviation from scale-invariance which 
stems from the time-dependence of the Hubble constant 
during the inflationary period). 

There are good heuristic reasons to expect a scale- 
invariant spectrum of fluctuations to emerge from infla- 
tion. Since de-Sitter space is time-translation-invariant, 
one should expect the amplitude of the density fluctu- 
ations SM/M to be independent of the scale (labelled 
by the comoving wavenumber n) if measured at the time 
when the corresponding wavelength crosses the Hubble 
radius In during the inflationary period. Since micro- 
physics cannot change the physical amplitude of the mass 
fluctuations while the wavelength is larger than In, one 
therefore expects SM/M to be independent of n when 
measured at the time if(n) when the scale re-enters 
the Hubble radius in the post-inflationary Friedmann- 
Robertson- Walker period: 

— {n,t{{n)) = const. , (1) 

which is the definition of a scale-invariant Harrison- 
Zel'dovich spectrum 

The time-translation-invariance is, however, broken in 
the current models of inflation. The calculations are done 
by picking an initial time ti (e.g. the beginning of the 
inflationary period), by choosing a specific state of the 
quantum fields at this time (e.g. the local Minkowski 
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vacuum state JTo) ] or the Bunch-Davies vacuum Jl2|), by 
evolving this state using the linearized equations of mo- 
tion, and by finally calculating the correlation functions 
and expectation values of interest. In this context, the 
emergence of a scale- invariant spectrum of fluctuations is 
seen to arise from a subtle cancellation of the wavenum- 
ber dependence in the initial state wave function and in 
the growth factor before Hubble radius crossing, and thus 
depends explicitly on the initial state chosen. States can 
be found |13|] which do not yield a scale-invariant spec- 
trum. Thus, it is clear that the prediction of a Harrison- 
Zcl'dovich spectrum is not completely generic in current 
models of inflation. 

There is, however, a much more serious potential prob- 
lem for the claim that current models of inflation based 
on weakly self-coupled scalar fields generically lead to a 
scale-invariant spectrum of fluctuations. Most of these 
models of inflation involve (see e.g. |l4| for a recent re- 
view) a period of inflation much longer than the 60 e- 
foldings of inflation required to solve the horizon and flat- 
ness problems of standard cosmology. Since wavelengths 
exponentially redshift during inflation, the physical wave- 
lengths of the modes which correspond to the present 
large-scale structure in the Universe were, in those mod- 
els, much smaller than the Planck length at the initial 
time t{. Thus, the usual computations of the spectrum 
of fluctuations involve extrapolating weakly self-coupled 
field theory coupled to classical gravity into a regime 
where these theories are known to break down. 

This problem is analogous to the Trans-Planckian 
problem for black hole physics (see Ref. for a re- 
cent overview). In black hole physics there is an arbi- 
trarily large blue shift when following modes of Hawking 
radiation at future infinity into the past, and the usual 
calculations of Hawking radiation |L6| seem suspect (see 
e.g. Ref. |l7j for a discussion of this point). 

In the case of the black hole problem, it was recently 
shown by Unruh |Q, Brout et al. Jigj] , Hambli and 
Burgess [^cj and by Corley and Jacobson that the 
prediction of a thermal Hawking spectrum of black hole 
radiation is insensitive to modifications of the physics at 
the ultraviolet end of the spectrum. In these works, the 
dispersion relation of the quantum fields was modified 
(in rather ad-hoc ways) at energies larger than some ul- 
traviolet scale kc , and it was found that the spectrum of 
radiation at future infinity at wavenumbers much smaller 
than kc is insensitive to the modifications considered. In 
this sense, Hawking radiation from black holes was shown 
to be an infrared effect. 

The obvious question is whether a similar conclusion 
will hold for the generation of fluctuations in inflation- 
ary cosmology. This is the question we will address in 
this paper. We will consider a free scalar field in an 
inflationary background [de Sitter phase of a Friedmann- 
Robertson- Walker cosmology with scale factor a(t)] . This 
scalar field can represent the scalar metric fluctuations, 
the gravitational wave mode, or a matter scalar field on 
the fixed background geometry - the case of most interest 



for cosmology corresponds to scalar metric fluctuations. 
We will modify the usual dispersion relation 

u 2 = k\ fc 2 ^' ( 2 ) 

where n and k are the comoving and physical wavenum- 
bers, respectively, for values of k larger than some cut- 
off scale fee j aud will calculate the predicted spectrum 
of fluctuations in the modified theory for well-motivated 
initial quantum states, states which in the unmodified 
theory coincide with the state usually chosen as the ini- 
tial state. The modified dispersion relations which we 
use are the same as the ones used by Unruh [[l8) and 
by Corley and Jacobson pi] . As preferred initial states 
we will use either the state which minimizes the energy 
density at the initial time U, following the approach of 
Brown and Dutton p5| ] , or a naive generalization of the 
local Minkowski vacuum. 

We find that in the case of Unruh's dispersion rela- 
tion, the spectrum of density fluctuations is unchanged 
in the minimum energy density initial state. However, in 
the case of the family of dispersion relations generaliz- 
ing the choice of Corley and Jacobson, the choice of the 
minimum energy density initial state leads to a spectrum 
of fluctuations which, depending on the specific mem- 
ber of the family of dispersion relations chosen, may be 
characterized by a tilt, by an exponential factor, and by 
superimposed oscillations. 

Our work indicates that the prediction of a scale- 
invariant spectrum in inflationary cosmology depends 
sensitively on hidden assumptions about super-Planck- 
scale physics. This has important implications for the 
attempts to unify fundamental physics and early Uni- 
verse cosmology. It is now a rather nontrivial question 
under which conditions a unified theory of all forces such 
as string or M-theory will lead to a scale-invariant spec- 
trum, assuming for the moment that it does indeed lead 
to a period of inflation. 

The outline of this paper is as follows. In Section II 
we demonstrate that the growth of linear density fluctu- 
ations, gravitational waves and linear scalar matter fluc- 
tuations can all be described in terms of the same frame- 
work: that of a free scalar field with a time-dependent 
mass. In Section III we introduce the two classes of mod- 
ified dispersion relations which will be used in the calcu- 
lations. The quantization of the scalar field in the time- 
dependent background and the construction of the min- 
imum energy density initial state are reviewed in Section 
IV. Section V contains our calculations for both classes 
of dispersion relations. Our results are summarized and 
discussed in the final section. 
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II. EQUIVALENCE BETWEEN COSMOLOGICAL 
PERTURBATIONS AND A FICTITIOUS 
SCALAR FIELD 

Without loss of generality, the line element for 
the spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLKW) background plus the perturbations can be writ- 
ten in the synchronous gauge according to [{24|]2E 



fluctuations of the inflaton are not coupled to fluctuations 
of the perturbed metric. Coupling occurs only as a result 
of the violation of the condition p = —p. 

Observable quantities can be computed when the ini- 
tial power spectra are known. These are defined in terms 
of the two-point correlation functions. For the Bardeen 
potential one has 



ds 2 =a 2 {r))\ -dr) 2 + 



h(rj,n)QSi. 



(0|$ B (J7,x)$ B (»7,x + r)|0) 
+ °° dn sin i 



+hi(rj, n)— ^ + h gw (r],n)Q i j 



-n 3 P $B (77,n), (6) 



dr'drf . (3) 

whereas for gravitational waves the correlator is given by 



In this equation, the dimensionless quantity n is the co- 
moving wavevector related to the physical wavevector k 
through the relation k = n/a{rf). r\ is the conformal 
time related to the cosmic time t by dt = a(rj)drj. The 
functions h and hi represent the scalar sector and Q{x l ) 
is the eigenfunction of the Laplace operator on the flat 
spacelike hypersurfaces. The function /i gw represents the 
gravitational waves and Qij{x l ) is the eigentensor of the 
Laplace operator. It is traceless and transverse, namely 



Qi — Qi 



0. It is convenient to introduce the back- 



ground quantity 7(77) defined by 7 = —H/H 2 , where 
a dot means differentiation with respect to cosmic time 
and H is the Hubble rate, H = a /a. We can also write 
7 = 1 — H' /H 2 , where H = a' /a and a prime denotes 
differentiation with respect to the conformal time. 

In the tensor sector, we define the quantity pt by 
h sv , = /it /«• Then, the equation of motion is given by 



/4 + 



a 
a 



"T = 0. 



(4) 



Since gravitational waves do not couple to matter, the 
last equation is valid for every type of matter. 

In the scalar sector, it is convenient to work with a 
residual gauge invariant variable ps defined by ps = 
[a/(W^/7)](/i' + Hjh) where we have supposed 7^0. 
The case 7 = must be treated separately (see be- 
low). The quantity ps is related to the gauge invari- 
ant Bardeen potential by $ B 3G) = [H-f/(2n 2 )}[p s /(a v ^f)}' 
where the subscript 'SG' means 'calculated in the syn- 
chronous gauge' p7| . Therefore, knowing the solution 
for /ig permits the calculation of the Bardeen variable. If 
matter is described by a scalar field (the inflaton), then 
one can show that p$ obeys the equation: 



// 

Ms 



Ms = 0. 



(5) 



The case 7 = corresponds to a scale factor a(t) oc e Ht , 
i.e. to the de Sitter manifold. Then, one can show that 
the exact solution to the perturbed Einstein equations 
is $b = 0: there are no density perturbations at all. 
This is because when the equation of state is p = —p, 



Jo n nr 

(7) 

where we have written = h gw Qij. We are specially 
interested in modes which are outside the horizon at the 
end of inflation, i.e. n/(aH) <C 1. For these modes, the 
power spectra do not depend on time and can be written 
as 



n 3 P$ B (n) = Agn 



n s — 1 



n 3 P h (n) = A T rf 



(8) 



Let us now consider power law inflation models where 
the scale factor is given by a{rf) = ZoM 1+/3 where /3 is a 
number such that [3 < —2 and Iq has the dimension of 
a length. The advantage of this class of models is that 
everything can be calculated exactly. In the case (3 = — 2 
which corresponds to exponential expansion, the length 
Iq is nothing but the Hubble radius, = a 2 /a' . The 
function 7 is a constant given by 7 = (/3+2)/(/3+l) which 
vanishes for (3 = —2. We see that Eq. (^) now reduces to 
Eq. (||). The spectral indices can be determined exactly 
and read 



n s = 2(3 + 5, 



tit 



2(3- 



(9) 



We have the relation ns — 1 = which is valid exactly 
only for power law inflation. 

Let us now consider a massless scalar field $(77, x) liv- 
ing in a FLRW spacetime. It is convenient to Fourier 
decompose the field and to introduce the quantity p de- 
fined according to $(r/,x) = [l/(2n) 3 ^ 2 ] J dn(/i/a)e m x . 
It is easy to show that the Klein- Gordon equation reduces 
to the following equation for p 



a 
a 



p = 0. 



(10) 



This equation is exactly the same as Eq. (|4[) and Eq. 
(||). Therefore, investigating the properties of cosmologi- 
cal perturbations is equivalent to investigating the prop- 
erties of a fictitious scalar field $(77, x). In particular, 
the calculation of the power spectrum of the scalar and 
tensor perturbations reduces to the computation of the 
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power spectrum of this fictitious scalar field. In the fol- 
lowing, we will restrict our considerations to this case, 
having in mind that, in fact, we will calculate the power 
spectra of cosmological perturbations. 

Let us make a last remark. Although it seems that we 
have considered only a limited class of models (i.e. power 
law inflation) , the previous analogy is in fact much more 
general. This is because the slow roll approximation, 
valid for a wide class of inflationary models, reduces to 
first order to power law inflation. 



III. TIME DEPENDENT DISPERSION 
RELATIONS 

In this section, we present the two classes of modified 
dispersion relations that will be used in this article. Let 
us return to the equation of motion (^|) . In this equation, 
the presence of the term n 2 is due to the differential op- 
erator 8 V didj in the Klein-Gordon equation. In Fourier 
space, this means that 



and neglect the term a" /a. Therefore, the solution can 
be expressed as: 



/in (77) = Bxe* 



(14) 



Finally, in Region III, the mode is outside the Hubble 
radius: A>(h and the solution (the growing mode) is 
given by: 



(15) 



where C is a n dependent constant. This constant has 
to be determined by performing the matching of /i and 
fi' at the times of transition between regions I and II 
and regions II and III, r/i and rj2 respectively. Then, the 
spectrum can be calculated and reads 



n 3 |C| 2 



(16) 



Let us now turn to the first example of a time dependent 
modified dispersion relation. 



to 



(11) 



The dispersion relation is therefore linear in the physi- 
cal wavenumber fc: lu = fc. A possible alteration of the 
high frequency behaviour of the Klein-Gordon equation 
can be obtained if we require the presence of a nonlinear 
function F(k) such that to = F(k) which, for physical 
wavenumbers smaller than a new characteristic scale fcc, 
i.e. k <C fcc, reduces to w » fc. This means that the 
n 2 term in the Klein-Gordon equation should now be re- 
placed with a time dependent n 2 s (rf) such that 



n 2 eS = a 2 ( V )F 2 (k) = a 2 (r,)F 2 [n/a(r,)]. 



(12) 



We see that, in terms of comoving wavenumbers, we ob- 
tain a time dependent dispersion relation. In what fol- 
lows, we will consider two explicit examples for the func- 
tion n c fj. Given the modified dispersion relation, Eq. 
(|l0|) can now be written as 



l cff 



a 
a 



fi = 0. 



(13) 



Let us analyze this equation in more detail. We can dis- 
tinguish three regimes. In Region I, the wavelength of a 
given mode, A(ry) = (2ir/n)a(ri), is much smaller than the 
characteristic length: A <C lc- The nonlinearities in the 
dispersion relation play an important role and the solu- 
tion of the equation of motion depends on the particular 
form of F(k). A crucial issue is that the mode no longer 
behaves as a free wave initially. As a consequence, the 
choice of initial conditions cannot be done in the usual 
way. In Region II, the wavelength of the mode is larger 
than the characteristic length but still smaller than the 
Hubble radius, lc <C A <C In- In this case, one can con- 
sider the dispersion relation to be linear, i.e. £l(rj) w 



A. Unruh's dispersion relation 

The dispersion relation used by Unruh in Ref. fL8|| , in 
the context of black holes physics, is: 



to 



F{k) = fc c tanh 



i/p 



fc 

fcc 



(17) 



where p is an arbitrary coefficient. For large values of the 
wave number, this becomes a constant fcc whereas for 
small values this is a linear law as expected. According 
to Eq. (|lj), in the context of cosmology, we take 



n cff (,) = ^Mtanh^ 



nlc V 
2ira(<o)J 



(18) 



where lc is the characteristic length corresponding to 
fcc- The argument of the hyperbolic tangent can also 
be rewritten as Ic/X{t)). This means that when A 3> lei 
n Q s(ji) tends to n. 



B. Generalized Corley/Jacobson dispersion relation 

The dispersion relation utilized by Corley and Jacob- 
son in Ref. 0] is given by the following expression 



uj 2 = F 2 (k) ^k 2 - w . (19) 
In this article, we consider a more general case and write 



= k 2 + k 2 J2b q \T- 



9 =1 



2<l 



(20) 



4 



where the b q are, a priori, arbitrary coefficients. Let 
us suppose that the previous sum only contains the last 
term. The physics depends on the sign of b m . If b m 
is negative, then u> vanishes for k = kc\b m \~ 2m . Beyond 
this point, the dispersion relation becomes complex. The 
Corley/Jacobson case corresponds to m = 1 and b\ = — 1 . 
In the context of cosmology, the previous ansatz gives rise 
to the following function n c s{rf) 



=i (2*0 



2q 



2q 



,29 



(21) 



Again, when A lc then the effective comoving 
wavenumber simply reduces to n. On the other hand, 
when A <C lc , one has 



2tt) 



2 m 



2m+2 



(22) 



The different dispersion relations used in this article 
are displayed in Fig. (|l|) together with the dispersion 
relation considered in Ref. [E2| denoted "KG" . 
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Jacobson/Corley b m >0 




Standard 








Unruh 1 




Jacobson/Corley b m <0 : 



k 



FIG. 1. Sketch of the different dispersion relations. 



IV. QUANTIZATION OF A MASSIVE SCALAR 
FIELD 

The aim of this section is to develop a Lagrangian and 
Hamiltonian formalism for the system described above. 
We will show that considering a time-dependent disper- 
sion relation is equivalent to giving a time-dependent 
mass to the fictitious scalar field. The main purpose of 
this section is to discuss the initial conditions. As already 
mentioned previously, when the wavelength of a mode is 



smaller than the critical length lc, the mode does not 
behave as a free wave because the dispersion relation in 
this region is no longer u> — k. As a consequence, it is no 
longer possible to impose the usual initial condition at 
77 = rji, i.e. fi — > e~ m ( r '~ T "V \/2~n. Another method must 
be used. Following Ref. |23|j . we will choose the state 
which initially minimizes the energy density of the field. 



A. Lagrangian and Hamiltonian formalisms 

We now study a massive fictitious scalar field $ whose 
action is given by 



S = 



,'2 



a z 



drj / dn 



(23) 



In this equation, the scalar field has been Fourier ex- 
panded according to 



$(77, x) 



1 



1 



(2tt) 3 / 2 a(r)) 



(24) 



and fJ-n(v) denotes the complex Fourier component of the 
field. We can easily check that the Lagrange equation of 
motion for the quantity /J. n (v) leads to Eq. (Jl3|). 

We are now in a position where we can pass to the 
Hamiltonian formalism. Our first move is to perform the 
following time-dependent transformation 



1 



N(n, rj) 



V>n(»7), 



(25) 



where N(n, rj) is a time-dependent factor which will be 
fixed below. Next, the action given in Eq. ( |23| ) expressed 
in terms of the new variable ipni 7 ]) takes the form 



S = 



~N 



N 2 \ N 



(26) 



We can now calculate the conjugate momentum to ipnirj). 
Its definition is p n = dC n / 'dip^irj) where C n is the La- 
grangian density (the bar indicates that one calculates 
the Lagrangian in Fourier space) which one can deduce 
from the previous equation. The conjugate momentum 
reads 

If a' \ N' 
Pn = ^(<--^J-^n. (27) 

The Hamiltonian can be determined using the following 
relation 
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Hn = Pn^n + Pn^'n ~ ^ 



(28) 



Inserting the expressions of the Lagrangian and of the 
the conjugate momentum in this definition, we obtain 



Hn 



i 



MO' 

aN 



N 2 



(29) 



«tl = 



The explicit quantization can now be carried out. We 
express the Fourier component tj} a and its conjugate mo- 
mentum p n in terms of creation and annihilation opera- 
tors, satisfying the usual commutation relation [c n 
6(n — r), according to 

ijj n = Vh\ r, 



Pi, 



Vhf 



(30) 



The Hamiltonian operator is obtained by plugging the 
previous expressions into Eq. (^9|) and requiring that 
l Ttu>/2' be present in each mode, which fixes the normal- 
ization factor N to be 



N 2 



Mv) , 



(31) 



where u> is the 'comoving frequency' defined by io(rj) = 
n e s- Although we use the same notation for convenience, 
this frequency should not be confused with the phys- 
ical frequency which appears in Eqns. ([l7j ) and (|l9| ) 
and which can obtained by multiplying the comoving fre- 
quency by a factor X/a. The Hamiltonian reads 



H 



dn 



ir- 



~2 



t t 



(32) 



This Hamiltonian has the usual structure. The first term 
is just a collection of harmonic oscillators whereas the 
second term represents the interaction between the back- 
ground and the perturbations. This term is responsi- 
ble for the phenomenon of particle creation, which is a 
squeezing effect. In a static spacetime, the pump func- 
tion (ay/uj)' I (ay/uj) vanishes and the interaction part of 
the Hamiltonian disappears. The field operator can be 
expressed as 



$(77, x) = 



1 



dn 



o(t?) (2tt)3/2 J v /2^) 



c„(77)e in 



(33) 



The time evolution of the creation and annihilation op- 
erators and therefore of the quantum scalar field is cal- 
culated by means of the Heisenberg equation: 



ih—Cn(r)) = [c n ,H]. 
arj 



(34) 



Using the form of the Hamiltonian derived previously, 
one gets the following equations of motion 



in— — = nui(r))c n + in i^c_ 



(35) 
(36) 



The solution of these equations is a Bogoliubov transfor- 
mation which can be written as 



dr\ 

ih—?- = -Hu)(ri)cl + ih K v ' c_ n . 
d?7 ay/uj 



(37) 
(38) 



where we have introduced two new functions u n (rj) and 
v n (r)). These functions satisfy \u n {rj)\ 2 — |^n(?7)| 2 = 1 in 
order for the commutation relation given to be preserved 
in time. Let us notice that u n and v n do not depend on 
the vector n but only on its modulus n. Inserting the 
previous equations in Eqns. (|35| ) and (|36|), one obtains 
the equation of motion for these two functions 



in— — = Jwvq)u n + in i=-v n , 

dl] dy/UJ 

. % dv n -AaVu)' » 
in— — = huj{r])v n + in 



(39) 
(40) 



The functions u n and v n can be re-expressed in terms of 
three other arbitrary functions r n (j]), 9 n (r]) and tp n (rj). 
Following this path would lead to the squeezed state for- 
malism. However, we will not need it in this article. 



B. Fixing the initial conditions 

The previous considerations permit to fix the initial 
value of the mode function /i n (?7i) an d its derivative 
MnC^i) f° r an y choice of function £1(t)), i.e. for any time 
dependent dispersion relation. 

It is straightforward to check that the function 



I 



N{n,rj) 



(u n + v*) 



1 



K + <), (4i) 



Mr, 



satisfies Eq. (|13J)Q From Eqns. (|37j) and we see 

that the initial conditions for the two function u n and 
v n are given by: u n (r) = = 1 and v n (r) = r/i) = 0. 
Therefore, the initial value of the mode function /1 can 
be written as: 



1 



1 



y/2u(rji) %/2n c ff 



(42) 



*It should be noticed that [i n is not exactly the mode func- 
tion introduced before. It is dimensionless (instead of dimen- 
sion s/Kc) and depends only on the modulus n. In the same 
manner, we now deal with a 'new' function ijj n = N/j, n . 
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Let us now turn to the determination of p'(r) — 771). It 
will be found by the requirement that the energy density 
is minimized. The stress energy tensor can be obtained 
from the action ( ^3|) with the help of the standard defini- 
tion. In terms of the Fourier components ip n , the energy 
density reads 



An 2 a 4 



-,'2 



-i>n1pn 



dn 

TV 2 



(aN)' 
aN 



(43) 



We now define the functions x(rj) and 2/(77) as the real 
and imaginary parts of the ratio ip' n /ip n = x + iy, respec- 
tively. Then, the initial energy density can be expressed 
in terms of X\ = x(n = rji), y\ = y{rj = rji) and the Wron- 
skian W(n) = p' n p^ — p„ p n which is a time independent 
quantity (as can be checked in calculating dW(n) / dr/ and 
using the equation of motion for p n ) 



47r 2 a 4 



'cff ' 



dn 





N' 2 

iV 2 " 



W(n) 
a'N r 



Vi 



„ (aN)' 



aN 



(44) 



where N and a are also evaluated at the initial time. No- 
tice that, while deriving the previous equation, we used 
the fact that the Wronskians of p n and ip n are related by 
a factor iV 2 . The 'vacuum' used in this article is defined 
as the state which initially minimizes the energy density. 
The variation of the previous expression with respect to 
Xi and y\ leads to 



S P : 



47r 2 a 4 



dn 



W(n) 
2i 



1 

yf 

N' 2 
TV 2 " 



Hi 



2^* 



(aN)' 



aN 



Sxi 



aN 



a'N' 



aN 



8y\ 



(45) 



Demanding that 5p = 0, one deduces the initial values of 
x and y 



a' . N' 
Xi = — (Vi) + -jrfW, yi = ±n e g. 
a N 



(46) 



These expressions can be simplified. Using the explicit 
form of the function N(n, rj), one can write 



N' _ J 
~~N 



(47) 



At the time r\ = rji, it is reasonable to consider A -C lc 
(otherwise, the whole problem studied here would be 
pointless). Then, for Unruh's dispersion relation, one 



finds N' /N w a' /2a and for the Corley/Jacobson disper- 
sion relation, one has N'/N ss —ma' /2a. In addition, 
a' /a is very small in the limit where the conformal time 
goes to —00 since a! /a(n\) = (1 + P)/\r]i\ and \rjy\ ^S> 1. 
Therefore, one gets that ip' n /ip n w iyi. On the other hand, 
we have p n = tj) n /N . Combining this formula with the 
previous one, one obtains p' n +N' /N = iyip n . Neglecting 
again the term N'/N, we finally arrive at 



p'(v = Vi) = ±« 



n c ff 
2 



(48) 



The initial conditions are now completely fixed and given 
by Eqns. © and ©. 

Let us also mention that it is possible to adopt another 
choice of initial conditions which corresponds to the 'in- 
stantaneous Minkowski vacuum' at 77 = rji, namely 



1 , jn 



(49) 



If the dispersion relation is standard, then £1 = and 
the mode is initially free: locally, it does not feel the 
curvature of space-time and behaves as it were flat. In 
this case, the two possible choices of initial conditions 
discussed above coincide. 

Two last comments are in order before ending this sec- 
tion. Let us first remark that the concept of an initial 
state which minimizes the energy density of the field 
could be problematic in a region where the dispersion 
relation becomes complex, as it is the case for the Cor- 
ley/Jacobson dispersion relation with b m < 0, since the 
energy needs not to be bounded from below in such a 
situation. We are not aware of any more obvious method 
than the one used here to deal with this case. 

Finally, although we have introduced two initial states, 
it should be clear that the minimizing energy state is 
the only physical vacuum state. The instantaneous 
Minkowski vacuum is considered here only to stress the 
fact that the choice of the initial conditions becomes 
more crucial than in the standard situation where one 
can show that a large class of initial states leads to the 
same spectrum p8[ (although, as already mentioned in 
the introduction, it is possible to find examples which do 
not belong to this class of initial states p3|). 



V. ANALYTICAL SOLUTIONS 

In this section, we calculate the spectrum of fluctua- 
tions for the two classes of dispersion relations introduced 
in Section III. We focus on a fixed comoving wavenum- 
ber 7i and proceed as follows. We solve the equation of 
motion in each of the three regions (defined in Section 
III) separately. The coefficients of the two fundamental 
solutions in Region I are fixed by the initial conditions 
discussed above. Then, we explicitly perform the match- 
ing of p and p! at the transitions between Region I and 
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Region II, which occurs at a time denoted by 771 , and 
between Region II and Region III, which occurs at time 
772, to obtain the coefficients of the two fundamental so- 
lutions in Region III, from which the spectrum can be 
calculated. The time r] 2 is when the mode crosses the 
Hubble radius, which is given by 



la 



\1 + P\ 



2+/3 



Thus, the condition luiv?) = -^(^2) boils down to 
, 2tt, 



\V2\ 



\l+p\ 



(50) 



(51) 



The geometry of space-time is illustrated in Fig. 




elenqth X = 2r, 



ia(7)) 



FIG. 2. Sketch of the evolution of a physical mode through- 
out the three regions denned in the text. The dashed region 
is the region where the dispersion relation is modified. 

We start this section with Unruh's dispersion relation. 



A. Unruh's case 



Note that, in fact, the form of this last equation is inde- 
pendent of the precise form (i.e. the hyperbolic tangent) 
of the dispersion relation in the regime lc A. It is 
just necessary to assume that F(k) goes to a constant. 
We see that the result depends in an essential way on 
the dimensionless parameter e = Ic/Iq- At this point, 
we have assumed nothing about the value of the ratio 
lc/lo- However, physically, it is clear that e <C 1. One 
would expect the cutoff length to be given by the Planck 
length (^c ~ ^pi), whereas ~ 10 5 Zpi if the spectrum of 
fluctuations is COBE normalized. In this case, we have 
e w 1CP 5 . In the following, we will use an expansion in 



terms of this parameter. The exact solution of Eq. ([55 
is 

f i 1 (r ] )^A 1 \f 1 r+A 2 \r ] r, (54) 
where the exponents x\ and x 2 are given by 



1 1 / 16tt 2 



(55) 



It is now time to fix the coefficients A\ and A2. T hey 
are completely determined by the initial conditions (|42| ) 
and (|48[). In the approximation where Iq ^> A, they are 
solutions of (note that we do not yet use the fact that e 
is small) 



Axfor +A 2 \ m \ x * = y^h\ 1/2 , 

A^hp- 1 +A 2 x 2 \r ]i p- 1 = T *\f^W-' 1 



(56) 
(57) 



The exact solution of this system of equations can be 
written as 



A x = \m\ 1/2 - Xl 

A 2 = h\ 1/2 - X2 



1 1 jef 2m 
1 — xx/x 2 2\ n\ ex 2 



1 



(58) 



y-(i±-). <*>. 

1 — x 2 /Xi 2 V 7T V exx 



The equation of motion for the mode function can be 
written as 



M " +< j^_ a 2 tanh 2/p 



( 1 



c 



a 
a 



/i = 0. (52) 



This equation can be solved exactly in Region I only if 
the scale factor is given by 0(77) = Iq/\t}\, i.e. in the 
case P — —2. Fortunately, this corresponds to the de 
Sitter space-time, the prototypical model of inflationary 
cosmology. Note that in this case Iq is the Hubble radius 
[see Eq. (p0|)]. In Region I, the hyperbolic tangent is 
approximatively one since ^ A initially. Therefore, 
Equation (|52|) reduces to 



\ v 2 



(i = 0. 



(53) 



It is at this point that we use the fact that e is small. To 
first order in a systematic expansion in this parameter 
we obtain 



x exp 



3/2 



1 2i7T 2lTT 

2-~ ± - 



2iir 



3/2 



• In 77; 



1 2i7T 2lTT 

2 + ~ ± - 



x exp [ In \rji 



(60) 



(61) 



We now pursue the calculation for both choices of the 
sign of the initial conditions. We introduce an index 'u' 
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for the upper choice and '1' for the lower choice. This 
leads to 



i 




16 








i( 


i) 




i) 







.3/2 



exp 



2in . 



■ In Vi 



1/2 



2z7T I 

exp ( In 1 771 



1/2 



- 1 

16 V 7T 



exp 

3/2 



In \r)i\ 

e 1 



exp [ In | ?/i 



(62) 
(63) 
(64) 
(65) 



Therefore, one has A% ^> A™, A\ 3> j4 2 and only one 
branch of the solution j54] ) survives. Then, the solution 
in Region I can be expressed as 



u,i, \ 1 {2 i7: , V \ 



(66) 



Let us now turn to Region II. As already mentioned 
above, in this region, the solution is given by 



inr\ 



(67) 



The coefficients B\ and B 2 are determined by the match- 
ing of this solution with the solution ( |66| ) at the time r\\ . 
Continuity of and // yields 



inB 2 e 



= m ( e|?7xj 
2 \ i 
mflie^" 1 - inB 2 e~ mr > 



1/2 



2l7T , 

exp [ =F m 



1/2 



2l7T , 

exp [ =p In 



The solution can be found easily and reads 

1 . .s 



Bi 
B 2 



2\/n 



2in 


In 


m 


+ 2iir 


e 


m 




2in 


In 


m 


2in 


e 


m 





(68) 
(69) 

(70) 
(71) 



As a consequence, the solution in Region II also contains 
only one branch. 

Finally, we must solve the mode equation in Region 
III. As already mentioned, the non-decaying mode is 



/ini = Ca(rj) . 



(72) 



The coefficient C is fixed by the matching of the mode 
function when the mode crosses the horizon at r\ 2 0- One 
gets 



c = mim)-;— = : — 



(73) 



Therefore, regardless of the choice of the sign of the initial 
conditions, we have |C| oc 1/n 3 / 2 and as a result 



n 3 P$ oc n° . 



(74) 



We see that, when (3 = —2, the final answer is not 
changed compared to what is obtained without the mod- 
ification of the dispersion relation, i.e. we get a scale 
invariant spectrum n.g = 1 [see Eq. (g)]. 

We now discuss different initial conditions. We adopt 
the 'instantaneaous Minkowski' initial conditions given 
by Eqns. (p)|). Of course, the form of the solution in 
Region I is still the same but, now, the coefficients A\ 
and A 2 are different. The exact expressions for theses 
coefficients can now be written as 



Ai = ± 
A 2 = ± 



1 1 — X\ 



2 x 2 — x\ 



2 x-i 



x 2 



it 
it 



IX 2 

ix\ 



(75) 
(76) 



In the limit when the parameter e is small, an expansion 
of the previous expressions leads to the following formulas 



A x 



1 



2V2n 



\m\- 2iw/e - 1/2 , a 2 



i 



2J2n 



V, 



|2i7r/e-l/2 



(77) 



The result does not depend on the choice of the sign of the 
initial conditions. We see also another crucial difference 
in comparison with the previous case, see Eq. ( |60| ) and 
(|6l|): this time, the coefficients are of the same order in 
e. Therefore, the solution in Region I is now given by a 
cosine instead of by a pure phase 



1 



1/2 



cos | — In 



(78) 



The solution in Region II is still given by plane waves. 
The matching at time r\\ permits the calculation of the 
coefficients B\ and B 2 . They read 



1 

2n 



B 2 = — 

2n 



7T 



27ri 



exp —inr\\ 



In 



exp i7i7yi 



27T?; , 

In 

e 



(79) 
(80) 



Again, there is an important difference in comparison 
with the previous case: both coefficients are now non 
vanishing. The mode function in Region II can be ex- 
pressed as 



^To be more precise, we should take the decaying mode in 
Region III into account and match both /i and p! at time 772. 
This only changes the result by an unimportant constant of 
order one. 



77 / 

— — - cos nn 



2tt 



77771 



■In 



(81) 



The function is proportional to 1/ti instead of 1/i/n. The 
determination of the constant C proceeds as previously 
and leads to the spectrum 
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n Pa oc n 1 1 



2?r 2tt , 
— + — In 

e e 



2tt 
n?7i 



(82) 



A few remarks are in order here. Firstly, the difference 
between ( |74| ) and (|8|) demonstrates that the final result 
does depend on the choice of the initial conditions. Sec- 
ondly, the spectral index is now modified and is ng = 
instead of ns — 1 previously. Thirdly, oscillations in the 
spectrum are present. If n\ and n 2 are two wave numbers 
such that the argument of the cosine differs by a factor 
2irp where p is an integer then one has n 2 jn\ = exp(pe). 
This means that unless p is comparable to e _1 , n\ and 
n 2 are almost equal. Therefore, the oscillations are very 
rapid. 



B. The Corley/Jacobson case 



With the dispersion relation 
tion becomes 



Of) , the equation of mo- 



H +fjt 



\V 



1+0 



2<l 



a 
a 



= 0. (83) 



This equation is valid for any scale factor of the form 
a(rj) = ZoM + ■ Unlike in Unruh's case, we do not need 
to specify (3 = — 2. 

We now need to discuss the form of the solution in Re- 
gion I. This crucially depends on the sign of the coefficient 
b m . In the regime we are interested in, i.e. lc ^ ^(v0> 
one can retain only the dominate term and the dispersion 
relation can be written as 

' , \ 2m 



A 



(84) 



This means that if b m is positive, the dispersion rela- 
tion remains real. If b m is negative the situation is more 
complicated. For very small value of \b m \, the dispersion 
relation can remain real even in the regime lc 3> A(?7i). 
However, it seems a bit artificial to fine-tune the value of 
\b m \ such that this actually happens. Without this fine- 
tuning the dispersion relation certainly becomes complex. 
This last property should not be considered as a surprise. 
Indeed there exist many situations in Physics where com- 
plex dispersion relations appear. This is for example the 
case in Hydrodynamics when one describes the damping 
of a sound wave in a fluid due to viscosity |2!| . Then, the 
dispersion relation is given by k — lo/c + iauJ 2 where a is 
a factor which depends on the viscosity coefficients. In 
Cosmology, other examples are Silk damping or damping 
of density perturbations due to neutrino decoupling Q . 
In this paper, we choose to analyze both cases and write 
b m = s\b m \ with s = ±1. Then, from Eqns. ( f42| ) and 
(Eq), the quantities [ii{r){) and A*j(?7i) take the form 



,-1/4 




1/2-6/2 



(85) 



M ^i)=±^/^M- 1/2+b/2 , (86) 



where we have defined b and 7 (not to be confused with 
the function 7 used in Section II) by the following ex- 
pressions 



b=l-m(l+j3), 7 = J™ e m n m+1 . (87) 



b{2ir) 



From the expressions ( J85| ) and (|8q), we deduce 

ti(m)/M = ±is 1 ' 2 b 1 \m\ h - x ■ (88) 

This ratio will out turn to be important in the calculation 
of the various coefficients determined by the matching 
procedure. To go further, we need to treat the cases 
s = ±1 separately. 

1. The case s = — 1, b m < 

In Region I, the equation of motion for the mode func- 
tion reduces to 



(27T) 



2m 



\v 



1+0 



n = o. 



(89) 



For a negative coefficient b m , the exact solution of Eq. 
( |S9| ) can be expressed in terms of modified Bessel func- 
tions as follows 



Mi(r?) = AM 1 ' 2 ^) + A 2 \r,\ 1 / 2 K v (z), 



(90) 



where v = 1/(26) and where the function z(rj) is defined 
by the following expression z(rj) = j\i]\ b . The coefficients 
Ai and A2 are determined by the initial conditions given 
in Eqns. (85) and (|86|). These coefficients should satisfy 
the system of equations 



(91) 



A x I u {zi) + A 2 K v (zi) = \rn\ 1/2 /xife), 

- AJy-iizi) + A 2 K^ 1 (z l ) - |7?i| \ 2 V ifa), (92) 

70 



where z\ denotes the value of z(r\) at time r\ = rji. The 
exact solution for A\ and A 2 can be expressed as 



A 1 = 7 | ?? i| b - 1 / 2 )U i(77 i )^_i(z i ) 

76 ^1(771)^-1(21) 

A 2 = - f \ri i \ b - 1/2 vi(m)r»-i(zi) 

, fal 1 -* Mm) iM) 



76 A*i(>7i) Iv-i{zi) 



(93) 



(94) 



In the derivation of the previous expressions, we used 
the exact equation: [I v K v -i + l v -\K v ){z) = l/z. Since, 
when Iq 3> A(r/i), the argument z\ is large we can now 
rewrite these equations using the asymptotic formulas for 
Bessel functions of large arguments Notice that it 

is necessary to go to the second order in the expansion 
of the modified Bessel functions. We obtain 
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At 



1/2 



i 1/2 m(m)h\ b/2 ~ 1/2 t 



i± i± 



2v- 1, 



2~ 



1-6 



-4, 



2tt 



1/2 



/ 1/2 Mi(^i)hi| &/2 " 1/2 ^ 



1T1T 



1 -2v, 



2l 



1-6 



(95) 



(96) 



For the sake of completeness, we pursue the calculation 
for both choices of the sign of the initial conditions. Let 
us again use an index "u" for the upper choice and "1" 
for the lower choice. We obtain: 



A\ 
A u 2 
A\ 
Al 



1 



1/2 



7 1/2 Mi(^i)k| fc/2 - 1/2 e- 2 



1/2 



-6/2-1/2 



1/2 



Mi(»?i)l»7i 



-6/2-1/2 



1v-\ 
2 7 V2 ' 



2tt 



1/2 



(97) 
(98) 
(99) 
(100) 



The exponential factor always determines the behaviour 
of the coefficients for any power of 1 771 1 . This implies A" « 
A\ w 0. We also see the following crucial effect: it turns 
out that for one choice of the sign of the derivative the 
first term in the squared bracket in Eqns. ([)||) and ( p^ ) 
cancels whereas for the other choice it is no longer the 
case. This has as a consequence that the dependence on 7 
is not the same. Since 7 depends on n, the n dependence 
of A 2 and A\ is not the same. We have A 2 oc r )~ 1 A l 2 . 

The second step of the calculation is to perform the 
matching of the solutions at the time rj = rj 1 . This will 
allow us to express the coefficients B\ and B 2 in terms of 
the coefficients A\ and A 2 . In Region II, the solution is 
given by plane waves. Therefore, the coefficients B\ and 
B 2 are now solutions of the equations 



B ie mm + B 2 e-' mrn - 



= A 1 \r h \ 1 / 2 I v (z 1 ) 
+A 2 \i 1i \ 1 ' 2 K u {z 1 ), 

B ie mr * - B 2 e~ mr * = -^AiM 6 - 1 / 2 /^!) 



(101) 



4 2| 



Ai 



Vi\ 1/2 Iv(zi) 



ijb 1^-1(^1) 
1 ^ ~T7 — 7~ 



+--f\m\ 1,2 K v { Zl ) 



(104) 



Much simpler (approximate) formulas can be obtained if 
one notices that the argument of the Bessel function is 
a big number Z \ = 7 |^i| f ' 3> 1, essentially because e is a 
small number in realistic cases. A very simple estimate 
allows us to quickly check the validity of this approxima- 
tion. We take m = 1, = 1, (3 = —2.2 which would 
correspond to a spectral index of rig = 0.6 for power law 
inflation and e = 10~ 5 as already discussed in the previ- 
ous section. We can then estimate z\ for n — An which 
corresponds to the mode which re-enters horizon today 
and which consequently mainly determines the value of 
the CMB quadrupole anisotropy. We find z\ ~ 4.7 x 10 4 . 
Therefore, we can again use the asymptotic behaviour of 
the Bessel function to simplify the previous equations. 
Putting all these ingredients together, we find \ 

B 1 « ^^y /2 | 77l |V2-6/ 2e -in, 1 -, 1 -i fj (105) 

B 2 a 4r(^V Vl^-fc/V^-^. (106) 

As a consequence, the solution in Region II can be writ- 
ten as 

Mv) = A 2 (j^j \ m \ 1/2 - b/2 e- Zl cos(nn - nr ll - . 

(107) 

The last step of the calculation is to perform the match- 
ing at i] — r\ 2 when the mode leaves the Hubble radius 
(boundary between Region II and Region III). As already 
mentioned, in Region III, the non-decaying solution is the 
super-Hubble function given by 



Mm fa) = Ca(r]) . 



(108) 



Repeating the same procedure as for Unruh's case, the 
spectra can easily be calculated and read 



n 3 P| CX „2/3+4 e 2( Zl - Zl ) cog 2 J ^ _ nm 



(109) 



+^A 2 \ Vl \ b -V 2 K^ 1 (z 1 ), (102) 72 3 P£ocn 2 P +2 - 2m e 2 ^- z ^ cos 2 (nr, 2 - n m -?-). (110) 
in \ 4/ 



where z\ is the value of the function z{r\) at r\ = rjx. The 
exact solution of this system of equations can be easily 
found and reads 



We see that the spectrum depends explicitly on the initial 
conditions chosen. We can check that the tilt is correct 
by noticing that n 3 P$ oc A\, 7 oc n m+1 and using the 



A2 



Ai 



1/2 1 A 



6-1 Iv-l 



>l) 



Vi\ 1/2 K„( Zl ) 



n I v {zi) 
, b _ x K v - X (zi) 



K v {zx) 



(103) 



'In order to be able to neglect the terms proportional to A\, 
we make use of the fact that \rji\ 3> 
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relation between and A\ already mentioned. From 
now on, we concentrate on the lower case which corre- 
sponds to an unmodified tilt and study the expression of 
the corresponding spectrum in more details (for conve- 
nience, we drop the subscript '1'). First, as mentioned 
above, we see that the power-law part is not modified in 
comparison with the usual case, i.e. the spectral index is 
still ns = 2/3 + 5. Secondly, there are oscillations in the 
spectrum since the argument of the cosine can be written 
as (considering for simplicity that \b m \ = 1) 



V(l+/3) 



(2+/3)/(l+/3) _ ^ 



(111) 



However, contrary to Unruh's case with Minkowski ini- 
tial conditions, no logarithmic dependence is present. 
Interestingly enough, for j3 — —2, the oscillations dis- 
appear. The most important part concerns the expo- 
nential factor. The factor z\ — Z\ is equal to Z\ — z\ = 
7M b ( 1 -MVk b ) ~ lW\ b = Zi since we have > |?7i|. 
The factor z\ can be re-written in such a way that the 
dependence on n is explicit 



\Z\bm\ fML il-m(l+/3)rn+l 

6(2tt)™ 



(112) 



The important factor in this expression is e m |?7i| 1 m ( 1 +' 3 ) 
since the others ones arc of order one. It can be re- written 



as 



. I l-m(l+/3) _ 



(113) 



and must be considered as large since \r)i\ 3> 1 and 
lc/a(Vi) ~ Ic/^iVi) ^ 1) a * least for wavenumbers not 
too different from 2ir. This means that the influence of 
the exponential factor is dominant and is responsible for 
a huge increase of the spectrum at large n. This is illus- 
trated if we write the spectrum for (3 = — 2 and m — 1 



n P$ oc e 



An' 



(114) 



where 4 » 1. Such a spectrum is almost certainly in 
contradiction with observations. 

We end this subsection with the calculation of the spec- 
trum in the case where the initial state is the minimum 
energy density state. We restart from the exact expres- 
sions for the coefficients A\ and A2, see Eqns. ( ^3|) and 
(£§). Using Eqns. @, we have 



\m 



1 1—6 



±- 



k 



< 1, 



(115) 



since, initially, Iq 3> A(?7i). As a consequence, we can 
derive a compact approximate expression for the coeffi- 
cients A\ and A2 



A\ « 0, A 2 



1 



(116) 



Notice that these formulas are valid for any choice of the 
sign of /ij(?7i). The rest of the calculation proceeds as 
above and leads to (|6 m | = 1) 



X COS 



2tt|1 + /3| 



2tt 



(117) 



The main difference in comparison with the spectrum of 
the previous section is the presence of a modified tilt. 
The spectral index is now given by ns = 2/3 + 5 + m. 



2. The case s = 1, b m > 

When the dispersion relation is real, the solution in 
Region I can be expressed in terms of usual Bessel func- 
tions 



Mi(ry) = Atlnl^Mz) + A 2 \ v \^ 2 J_ v (z), 



(118) 



where v and z(r\) have already been defined previously. 
The coefficients A\ and A 2 are now solutions of the fol- 
lowing system of equations 

A l J v {z i ) + A 2 J_„{zi) = /ii(?7i)|r7ir 1/2 , 



AiJu(zi) + A 2 J^(zi) = 

70 



(119) 
(120) 



Using the relation expressing the Wronskian \3- v J v -\ + 
J-v+iJv]{z) = 2 sin[7r/(26)]/(7rz), and performing some 
straightforward algebraic manipulations, exact expres- 
sions can be easily found. They read 

7T7 



A x = 



2 sm(niy) 



\V 



16-1/2 



1 



Iml 1 -" /4fa0 j_ v { Zi ) 



Ao 



1 



76 /J>i(r]i) Ji-v(zi) 

ml"- 1 ' 



2 sm(7rz/) 



(121) 



(122) 



76 Ml(?7i) Jv-l(Zi)_ 

These expressions are not valid if v = 1 /(2b) is an integer 
and this particular case must be treated separately. In 
this article, we assume that this does not happen. Since 
the solution in Region II is still given by plane waves, the 
derivation of exact expressions for the coefficients B\ and 
B 2 proceeds as before. Explicit matching of the mode 
function and of its derivative leads to 



n J v (zi) 



+^\m\ 1/2 J-Az,) 



n J-v\z\) 



(123) 



e - mm B 2 



m\ 1/2 Mzi) 



1 



n Jv{Zi) 



\vi\ 1/2 J-Azi) 



n J-u(z\) 



(124) 
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Having all the relevant exact expressions at our disposal 
we can now start to do some approximations based on 
the fact that z\ is a big number. For convenience, we 
introduce two new definitions [not to be confused with 
the functions x(rj) and y(rj) introduced in section IV-B] 

x(r)) = z{rj) + y ^ = z ^> ~ T ~ 4' ^ ' 

Then, using the expressions of the Bessel functions for 
large arguments we find 



]b J2-y,J^± iXi (12g) 
2 J " 1 sin(TTiy) v ; 



A^^^h^-^J^e^, (127) 
\ 2 / sm(7r^) 

where a;j = 37(771) and y; = 2/(771) - The correct matching 
time is M = [ne/(2^)]^ 1+ ^ /[2m(1+/3)1 , see Ref. @, 
and is equal to the time at which A = lc if &m = 1- bi 
the following, for simplicity, we consider b m = 1. The 
coefficients B\ and B2 can be expressed as 

B 1 « (2^) 1/2 |^i 1 1 / 2 - fo /2 e — n cos yi 
— Ij4i sin j/i + cos xi — 1A2 sin zi I , (128) 



B 



1/2 



j^j 1/2-6/2^ f^ Ai 



cos yi 



\27r7 J 

+iA\ sinyi + ^2 cosxi + 1A2 sinxi I , (129) 



where x\ = x{rj\) and yi — y{rji). Our next move is to 
repla ce t he expressions of A\ and A2, see Eqns. (126) 
and (127), in the previous formula. This leads to 



B 1 = Ti 



inr\\ 



2 sin(-7w) 



1/2-6/2 



e ±iXi I cosyi 



-i sin yi — e Tl7Tl ' cos x± + «e =R7r " sin £1 



B 2 =Ti 



2 sin(7w) 



1/2-6/2 



• T 1TTU • T i~KU • 

-j sini/i - cosxi — ie + smxi 



(130) 



(131) 



Then, the mode function at time rj = i]2 (which is the 
relevant quantity for the determination of the constant 
C) can be expressed as 



2 sin(-7w) 



1/2-6/2 



(132) 



From this equation, the expression of the spectrum can 
be easily established and reads 



n P$ cx n 



2/3+4 



(133) 



Let us analyze this spectrum in more detail. The first 
remark is that the tilt is unchanged and that the spectral 
index is given by the usual expression ng = 2/3 + 5. The 
second remark is that the exponential dependence has 
disappeared. This is due to the fact that, for s = 1, this 
factor becomes a pure phase. We recover the usual result 
as pointed out in Ref. j3j| . 

Let us finally turn to the case where the initial con- 
ditions are those which correspond to the instantaneous 
Minkowski vacuum. Restarting from the exac t ex pres- 
sions for the coefficients A\ and A2, see Eqns. (121) and 
(|122|), and using Eqns. (E9I), we find 



(?) M^^sin^, (134) 
\ 2 / sm(7ri') 

ki\ 1/2 \ ,b/2-i/2 wO/O . n „ K , 

~7T J h\' ' — f — rsmyi. (135) 
2 / sm(7r^) 



At 
A 2 



Inserting these equations into the exact formulas giving 
the coefficients B\ and B 2 , we obtain 



Bi 



W(?7i)e" 



-6/2+1/2 



2 sin(7w) 

x sin Xi cos yi — 1 sin Xi sin j/i — sin j/j cos x\ 



+1 smy; smxi , 



(136) 



B 2 



-6/2+1/2 



2 sin(7ri/) 

x ( sin Xi cos yi + i sin Xi sin j/i — sin y^ cos 



-2 smy; smii 



(137) 



We are now in a position where we can write the expres- 
sion of the mode function at time r\ — 772 • It reads 



Ad(^i) 


m 


-6/2+1/2 r 


2 sin(7w) 


Vi 





B(n)e in(v2 - m) 



+B*{n)e- mi - ri2 -^ , 
where the function B(n) is defined by 



(138) 



B(n) = sin x t cos y\ — i sin Xi sin j/! — sin y\ cos a; 1 

+z sinyj sina;i. (139) 

Then, one can write B(n) as B(n) = \B\e 1 ^ and define 
B(n) as S(n) = |-B(n)| cos(nrj2 — ni]i +-0). It follows that 
the spectrum can be written as 



n^ocn^+H^n)! 2 . 



(140) 
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The spectral index is given by rig = 2/3 + 5 + m, i.e. it 
differs from the standard one but is equal to the spec- 
tral index obtained in the case s = — 1 for instantaneous 
Minkowski initial conditions. The factor |£?(n)| 2 is of or- 
der one and produces a complicated oscillatory pattern. 

In conclusion, the resulting spectrum in the case of 
the Corley/ Jacobson dispersion relation is very different 
from the usual spectrum calculated using an unmodified 
dispersion relation, and different from what is obtained 
using Unruh's relation, even for initial conditions which 
minimize the energy. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied the dependence of the predictions of 
inflationary cosmology for the spectrum of fluctuations 
on hidden assumptions about super-Planck-scale physics. 
The motivation for our work is that in most current mod- 
els of inflation, the period of exponential expansion lasts 
so long that at the beginning of inflation, scales of cos- 
mological interest today had a physical wavelength much 
smaller than the Planck length, and the theories used to 
compute the spectrum of fluctuations are known to break 
down on these scales. 

We studied the problem by replacing the dispersion re- 
lation of a free field theory which is used to compute the 
spectrum in the standard approaches by a modified dis- 
persion relation, the modifications only being important 
on length scales smaller than a cutoff length Iq (which 
we expect to be given by the Planck length) . We consid- 
ered two classes of dispersion relations, based on the ones 
considered by Unruh jl8j (Class A) and by Corley and 
Jacobson |gl| (Class B), respectively, in their studies of 
the trans-Planckian problem of black hole physics. Ad- 
mittedly, modifying the physics in this way is a very ad 
hoc way of taking into account possible effects of super- 
Planck-scale physics, chosen for mathematical simplicity. 
We do not want to introduce mode-mode coupling in or- 
der to keep the computations simple. However, in order 
to demonstrate that there is a possible problem for the 
robustness of the usual predictions of inflation it is suffi- 
cient to construct one example of a modified theory which 
leads to different predictions. 

For a non-standard dispersion relation the choice of 
initial state becomes more difficult. We considered two 
choices, both of which coincide with the usual initial state 
in the case of the standard dispersion relation. The first 
(and better motivated) choice is the state which mini- 
mizes the energy density, the second choice is the naive 
generalization of the 'local Minkowski vacuum'. 

We have shown that for Class A dispersion relations 
the usual predictions of inflationary cosmology are re- 
covered (in the case of exponential inflation) if the initial 
state minimizes the energy density. In particular, the 
spectrum of fluctuations is scale-invariant. If the initial 
state is chosen to be the 'local Minkowski vacuum', then 



the resulting spectrum has a tilt and superimposed oscil- 
lations. 

In contrast, for Class B dispersion relations and an 
initial state which minimizes the energy density, the re- 
sulting spectrum of fluctuations is in general not scale- 
invariant. The precise nature of the spectrum depends 
sensitively on whether the dispersion relation turns com- 
plex or remains real. In the complex case, the spectrum 
is characterized by an exponential factor (more power in 
the blue, i.e. ns > 1), a tilt (compared to the "standard" 
predictions) which depends on the precise initial condi- 
tions, and superimposed oscillations. The exponent, the 
tilt, and the precise oscillatory pattern depend on the 
specific member of the class of dispersion relations cho- 
sen. For a spectrum which remains real, the usual result 
is unchanged. 

The reason why for Class A dispersion relations the 
usual predictions of inflation are maintained is that the 
time evolution during the period when the mode wave- 
length is smaller than the cutoff scale is adiabatic. This 
emerges from our calculations, but an intuitive way of 
understanding the result is that at all times the effec- 
tive frequency of the mode is larger than the Hubble rate 
and the initial vacuum state therefore adjusts itself adi- 
abatically to track the instantaneous vacuum state, thus 
leading to the same state at time rji as in the theory with 
unmodified dispersion relation For Class B dispersion 
relations, in contrast, the dispersion relation varies too 
quickly as a function of time while the scale is smaller 
than the critical length Iq and hence the evolution is not 
adiabatic . 

Let us now be more quantitative about the previous 
discussion. In Region I, Eq. ( |l3| ) can be written as 



(141) 



where u> p h ys is the physical frequency defined by w p h ys = 
(l/a)v/n 2 + a 2 VL 2 {n,r 1 )/l 2 c . The latter can be considered 
as constant as long as its characteristic time scale of evo- 
lution is small compared to the Hubble time, i.e. as long 
as we have adiabaticity. Therefore, let us define an "adi- 
abaticity coefficient" a according to 



a(n,r}) 



H 



drj 



(142) 



where we recall that H = a' /a . When a> 1, adiabatic- 
ity is satisfied and Eq. (141) reduces to the equation 
of motion in the Unruh's case, Eq. (|53|). In this situa- 
tion, we know that the final spectrum is unmodified since 
there is an exact cancellation of the n-dependence in the 
minimizing energy state and in the growth factor before 
Hubble radius crossing which results in the usual spec- 
trum. The previous argument shows that an unmodified 



^We thank Bill Unruh for making this point to us. 
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spectrum is expected when a 1 in the region where 
the dispersion relation is modified. Let us also note in 
passing that for the standard case, a = 1, since the time 
scale of evolution of w p hy S (^j v) an d °f the Hubble rate is 
the same. 

We have calculated the adiabaticity coefficient for the 
different cases treated in this article. The result is dis- 
played in Fig. (||). When r\ goes to — oo, we have Iq ^ A 
whereas Iq <C A when 7/ goes to zero. 

We see that there exists a clear difference between Un- 
ruh's case and the Corley/Jacobson cases. In the Unruh's 
case, a goes to infinity when lc » A and adiabaticity is 
preserved. When b m < 0, the adiabaticity coefficient 
reaches zero at the time when w = . Then, adiabatic- 
ity is progressively re-established. The coefficient a goes 
to infinity and there is a divergence when duj/dn = . 
In the regime when Iq <C A, a goes to one as it should 
since the various dispersion relations all become similar 
to the standard one. The previous considerations ex- 
plain why the final spectrum can be modified in the Cor- 
ley/Jacobson case with a complex dispersion relation but 
not in Unruh's case. 
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FIG. 3. Sketch of the time evolution of the adiabaticity 
coefficient a for the dispersion relations. 

We conclude that it is possible that in models of infla- 
tion based consistently on a unified theory at the Planck 
scale the predictions for fluctuations will not coincide 
with the usual predictions from our current inflationary 
Universe models. Generalizing from our results, a crucial 
issue appears to be whether the evolution of the quantum 
states corresponding to the fluctuations will be adiabatic 
on length scales smaller than the Planck scale. 

Our results point to the possibility that the interac- 
tion between fundamental physics and cosmology may be 
much richer than hitherto assumed. It is not only a ques- 



tion of if and how fundamental physics leads to inflation; 
a much richer question is what the specific predictions 
of the fundamental model of inflation will be, assuming 
for the sake of argument that such a fundamental model 
exists. 

It is not surprising that super-Planck-scale physics may 
modify the usual predictions of inflation. One model of 
the early Universe motivated by string theory, the Pre- 
Big-Bang Cosmology |Q based on dilaton gravity, leads 
to a super-exponential period of early evolution in which 
the Hubble constant is increasing, and where the pre- 
dicted spectrum of scalar metric fluctuations is not scale- 
invariant pq ]. It would be interesting to analyze the 
predictions of other models of inflation based on string 
theory, taking into account the evolution on string scales. 
One toy model in which this question could be analyzed is 
the nonsingular Universe |3q] based on higher derivative 
terms in the gravitational action. 

In the context of the models studied here, it would 
be interesting to explore whether the minimum energy 
density initial state is an attractor in a similar sense that 
the local Minkowski vacuum is in standard inflationary 
cosmology p8[ . 

Note that models of inflation based on a strongly inter- 
acting theory (such as the model analyzed in p7[ ) do not 
suffer from the Trans-Planckian problem discussed in this 
paper. In strongly interacting theories, perturbations are 
generated at all times at a fixed physical scale, and a 
scale-invariant spectrum results based on the heuristic 
arguments mentioned in the Introduction. In such theo- 
ries, however, the presence of strong interactions makes 
it hard to calculate the amplitude of the resulting spec- 
trum. 
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